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By H. Schlichting 


1. INTRODUCTION 


The flow about a body of revolution rotating about its axis and 
simultaneously subjected to an airstream in the direction of the axis 
of rotation is of importance for the ballistics of projectiles with 
spin. In jet engines of all kinds , too, an important role is played 
by the flow phenomena on a body which is situated in a flow and which 
at the same time performs a rotary motion. Investigations of 
C. Wieselsberger-*- regarding the air drag of slender bodies of revolution 
which rotate about their axis and are at the same time subjected to a 
flow in the direction of the axis of rotation showed a considerable 
increase of the drag with the ratio of the circumferential velocity to 
the free-stream velocity - increasing more and more, the slenderer the 
body. Similar results were obtained by S. Luthander and A. Rydberg^ 
in tests on rotating spheres which are subjected to a flow in the direc- 
tion of the axis of rotation. These authors observed, in particular, a 
considerable shifting of the critical Reynolds number of the sphere 
dependent on the ratio of the circumferential velocity to the free-stream 
velocity. The physical reason for these phenomena may be found in the 
processes in the friction layer where, due to the rotary motion, the 
fluid corotates in the neighborhood of the wall and, consequently, is 
subjected to the influence of a strong centrifugal force. It is clear 
that the process of separation and also the transition from laminar to 
turbulent conditions are strongly affected thereby, and that, therefore, 
the rotary motion must exert a strong influence on the drag of the body. 


*"Die laminare Stromung um einen axial angestromten rotierenden 
Drehkorper . " Ingenieur-Archiv, vol. XXI, no. 4, 1953, PP- 227-244.- An 
abstract from this report was read on the VIII International Mechanics 
Congress in Istanbul on August 27, 1952 . 

•*€. Wieselsberger, Phys. Z. 28 , 1927 , p. 84. 

O 

S. Luthander, A. Rydberg, Phys. Z. 36 , 1935, P- 552. 
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In the flow processes in the corotating layer of the fluid, one deals 
with complicated three-dimensional boundary-layer flows which so far ha ve 
been little investigated, experimentally as well as theoretically. Th. v. 
Karman^ treated at an early date the special case of a disk rotating in a 
stationary liquid, for laminar and turbulent flow, as a boundary-layer 
problem, according to an approximation method. Later, W. G. Cochran** also 
solved this problem for the laminar case as an exact solution of the 
Navier -Stokes equations. A generalization of this case, namely the flow 
about a rotating disk in a flow approaching in the direction of the axis 
of rotation, for laminar flow, has been treated recently by H. Schlichting 
and E. Truckenbrodt5 . The result most important for practical purposes 
are the formulas for the torque of the rotating disk; it is highly depend- 
ent on the ratio of the circumferential velocity to the free -stream veloc- 
ity of the disk. 

For the general case of a rotating body simultaneously subjected to 
a flow, J. M. Burgers^ gave a few general formulations . We have set our- 
selves the problem of calculating the laminar flow on a body of revolution 
in an axial flow which simultaneously rotates about its axis7. The prob- 
lem mentioned above, the flow about a rotating disk in a flow, which we 
solved some time ago, represents the first step in the calculation of the 
flow on the rotating body of revolution in a flow insofar as, in the case 
of a round nose, a small region about the front stagnation point of the 
body of revolution may be replaced by its tangential plane. 

In our problem regarding the rotating body of revolution in a flow, 
for laminar flow, one of the limiting cases is known: that of the body 

which is in an axial approach flow but does not rotate. The solution of 

^Th. v. Karman, Z. angew. Math. Mech. 1, 1921, p. 255- 

k 

W. G. Cochran, Proc. Cambridge Philos. Soc. 50, 193**, P- 365. 

H. Schlichting, E. Truckenbrodt, Z. angew. Math. Mech. 32, 1952, 
p. 97; abstract in Journal Aeron. Sciences 18, 1951, p. 638. 

^ J . M. Burgers, Kon. Akad. van Wetenschappen, Amsterdam *+5, 19**1.> 
p. 13. 

7 It is pointed out that the turbulent case, for the rotating disk in 
a flow as well as for the rotating body of revolution in a flow, mean- 
while has been solved, in continuation of the present investigations, by 
E. Truckenbrodt. Publication will take place later.- E. Truckenbrodt, 

"Die Stromung an einer angeblasenen rotierenden Scheibe bei turbulenter 
Stromung," will be published in Z. angew. Math. Mech.- E. Truckenbrodt 
"Ein Quadraturverfahren zur Berechnung der Reibungsschicht an axial 
angestromten rotierenden Drehkorpern." Report 52/20 of the Institut fur 
Stromungsmechanik der T. H. Braunschweig, 1952. 
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this case was given by S. Tomotika®, by means of transfer of the well 
known approximation method of K. Pohlhausen9 to the rotationally symmet- 
rical case. The other limiting case, namely the flow in the neighborhood 
of a body which rotates but .is not subjected to a flow is known only for 
the rotating circular cylinder^- 0 , aside from the rotating disk. In the 
case of the cylinder one deals with a distribution of the circumferential 
velocity according to the law v = olfi/v where R signifies the cylinder 
radius, r the distance from the center, and to the angular velocity of 
the rotation. The velocity distribution as it is produced here by the 
friction effect is therefore the same as in the neighborhood of a poten- 
tial vortex. In contrast to the first limiting case (nonrotating body 
subjected to a flow), the flow in the case of slender bodies which rotate 
about their longitudinal axis in a stationary fluid does not have 
"boundary-layer character," that is, the friction effect is not limited 
to a thin layer in the proximity of the wall but takes effect in the 
entire environment of the rotating body. 

Very recently, L. Howarth-^- also made an attempt at solution for a 
sphere rotating in a stationary fluid. This flow is of such a type that 
in the friction layer the fluid is transported by the centrifugal forces 
from the poles to the equator, and in the equator plane flows off toward 
the outside. 

When we treat, in what follows, the general case of the rotating body 
of revolution in a flow according to the calculation methods of Rrandtl's 
boundary-layer theory, we must keep in mind that this solution cannot con- 
tain the limiting case of the body of revolution which only rotates but 
is not subjected to a flow. However, this is no essential limitation 
since this case is not of particular importance for practical purposes. 


The dominant dimensionless quantity for our problem is the ratio 


Circumferential velocity V m R^n 
Free-stream velocity Uqq 


where R m is to denote the radius of the maximum cross section of the 
body of revolution. The calculations must aim at determining for a 
prescribed body of revolution the torque, the drag, and beyond that, 
the entire boundary-layer variation as a function of V m j U OT . The 

®S. Tomotika, "Laminar Boundary Layer on the Surface of a Sphere 
in a Uniform Stream." ARC Rep. 1678 , 1935. 

9k. Pohlhausen, Z. angew. Math. Mech. 1, 1921, p. 253. 

10 H. Schlichting, "Grenzschicht-Theorie, " p. 63 . Karlsruhe 1951 . 

1 ' L L. Howarth, Philos. Mag. VII Ser. 42, 1951, p. 1,308. 
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particular case ’V m /^ oa = 0 is already known from the boundary-layer 

theory established so far. Considering what has been said above, we 
must not expect our solution to be valid for arbitrarily large 

The upper limit of the value of V m / U^, for which our calculations 

holds true, still remains to be determined. Presumably, it will lie 
considerably above V m ju m = 1 . 


2. THE FUNDAMENTAL EQUATIONS 


We take the coordinate system indicated in figure 1 as a basis for 
the calculation of the flow. Let (x,y,z) be a rectangular curvilinear 
fixed coordinate system. Let the x-axis be measured along a meridional 
section, and the y-axis along a circular cross section so that the 
xy -plane is the tangential plane. The z-axis is at right angles to the 
tangential plane. Let u, v, w be the velocity components in the direc- 
tion of these three coordinate axes. Furthermore, let R(x) be the 
radius of the circular cross section, co the angular velocity of rota- 
tion, U(x) the potential -theoretical velocity distribution, and 
v = p/p the kinematic viscosity. 

The equations of motion simplified according to the calculation 
methods of the boundary-layer theory are for this coordinate system 


^ + u dR + dw _ o (continuity) (l) 

dx R dx dz 


(momentum, meridional) ( 2 ) 


(momentum, azimuthal) ( 3 ) 


du v 2 dR chi dU d2u 

u + w — = U — + v — — 

dx R dx dz dx dz 2 


dx R dx dz dz 2 
The boundary conditions are 


z = 0: u=0, v = Vq = Rco, w = 0; 


u = U(x), v = 0 (ij-) 



NACA TM 114-15 


5 


A solution of this system of differential equations for an arbitrarily 
prescribed body shape R(x) with the pertaining potential theoretical 
velocity distribution U(x) leads to insurmountable mathematical diffi- 
culties.. We use therefore the more convenient approximation method which 
makes use of the momentum' theorem. We obtain the two momentum equations 
for the meridional and the azimuthal direction by integration of the 
corresponding equations of motion over z from the wall z = 0 to a 
distance z = h > 5 which lies outside of the friction layer. 

For the meridional direction there results by integration of (2) 
over z, with consideration of the continuity equation (1) and after 
introduction of the wall-shear stress for the x-direction 



the momentum equation for the meridional direction 


( 5 ) 


♦ V) ♦ I fpx ♦ *0%) - If (6) 

Therein, as is well known. 



is the displacement thickness whereas 



may be denoted as momentum-loss thicknesses for the x- or y-direction. 
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In an analogous manner there results for the azimuthal direction by- 
integration of ( 3 ) over z with consideration of the continuity equa- 
tion (l) and introduction of the wall-shear stress for the y-direction 



as the momentum theorem for the circumferential direction 


( 10 ) 


Therein 


d 

dx 


(ur 3 V ) , .s2lm 


xy 


Ja U 




0 u v 0 


(ID 


( 12 ) 


has been introduced as the "momentum loss thickness due to spin. 


3. APPROXIMATION METHODS 
(a) The Velocity Distributions 


According to the approximation method of the boundary-layer theory 
as given first by Th. v. Karman and K. Pohlhausen, the momentum equa- 
tions ( 6 ) and ( 11 ) are satisfied by setting up suitable formulations for 
the velocity distributions u and v which satisfy the most important 
boundary conditions. For the present case, two parameters may still be 
left undetermined in these equations for the determination of which the 
two momentum equations are then available. As expressions for the veloc- 
ity distribution, polynomials in the distance from the wall have proved 
to be suitable, with the property that the boundary layer joins at a 
finite wall distance z = 5 the frictionless outer flow. The boundary- 
layer thickness may be different for the meridional and the azimut ha l 
velocity component. Let these boundary-layer thicknesses be 6 ^- and 

S y , respectively; we introduce the dimensionless wall distances formed 
with them 


— = t and — = t ' 

5x 5y 


( 13 ) 
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For the velocity distributions u and v, we select polynomials of 
the fourth degree in t and t ' , respectively. These contain five coeffi- 
cients each so that , for determination of these coefficients, we can 
satisfy five boundary conditions each for u and v. We choose the 
following 10 boundary conditions : 


t = 0: u = 0 


t = 1 : u = U 


v a 2 u = _udU _ Vl M 

8z 2 dx R dx 

> 


= 0 & = 0 

dz 8z 2 


(l 4 a,b,c,d,e) 


t 1 = 0: v = Vq = Rfo 


SlS = o 

3 z 2 


t* = 1: v = 0 


^ 1=0 

8z 



(I5a,b,c,d,e) 


The boundary conditions (l 4 a, b, c) and (15a, b, c) result immediately 
from the fundamental equations (2) and (3) with ( 4 ) for z = 0 and 
z = 5 X or by. The remaining boundary conditions provide a gentle 

transition of the boundary layer into the outer flow. Taking these 
boundary conditions into consideration, one obtains the following poly- 
nomials as expressions for the velocity distributions 


— = 2 t — 2 t 3 + t^ + K— ft - 3 t 2 + 3 t 3 _ t^ 

u 6 ' 


(16) 


JL = 1 . 2t' + 2 t ’3 - 
v 0 


(17) 


Therein 


HU + 

(si 

.2 I 

\ U cffi 

dx 


/ R dx 


K = 


(18) 
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is a form parameter of the u-velocity profile , which is analogous to the 
form parameter A of the Pohlhausen method^- 2 . The velocity distributions 
u/u and v/vq are represented in figure 2. 


Let the point of separation be given by the beginning of the return 
flow of the meridional velocity component u(z) in the proximity of the 
wall 



This yields 


K = -12 (separation) 


(19) 


The expression for the u-component is the same as in the Pohlhausen 
method for the plane and rotationally symmetrical case. This guarantees 
that our solution in the case without rotation, to = 0 , will be trans- 
formed into the solution of S. Tomotika and F. W. Scholkemeyerl3 for the 
nonrotating body of revolution. Introduction of the expressions (l 6 ) 
and ( 17 ) into the momentum equations ( 6 ) and (ll) yields two differential 
equations for the still unknown boundary-layer thicknesses 5 x (x) and 

5 (x) or the quantities derived from them. 

o 


(b) The Momentum Equation for the Circumferential Direction 

We present first the further calculation for the momentum equation 
of the circumferential direction. With 



P 



-2— Rco 

Sy 


( 20 ) 


there results from (ll), after division by to. 



( 21 ) 


■^Cf. H. Schlichting, "Grenzschicht-Theorie, " p. 193* 

^F. W. Scholkemeyer, "Die laminare Reibungsschicht an rotations- 
symmetrischen Korpern." Dissertation Braunschweig 19^3, Cf. H. 
Schlichting, Grenzschicht-Theorie, p. 204. 
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With introduction of the further parameters 


- «=zie»s»o*:. 


as well as 



and 




and 


( 22 ) 


(23) 


CT = 


e (M + 5U 
\dx R d x/ 


( 24 ) 


one obtains from (21) the following differential equation for 0(x) 


d0 = G(K.A) 
dx U 

Therein is 

G(K,A) = 4g 0 - 2 a 

a universal function of the two parameters K and A. 

This function has been determined already by W. Dienemann-*-^ in the 
calculation of the temperature boundary layer on a cylinder (two- 

dimensional problem).^ p or the temperature distribution in the boundary 
layer there we chose the same polynomial of the fourth degree as we did 
for the azimuthal velocity distribution according to ( 17 ). According to 
(12) we have 


(23) 

(26) 


SO 


_ *»y __ 


Sy 


Because of 


u 


idf 


Jo u V Q 


t = 


t'A 

5 X 5 y 6 X 


Dienemann, "Berechnung des Warmeuberganges an laminar ums tromten 
Korpern mit konstanter und ortsveranderlicher Wandtemperatur . " Disserta- 
tion Braunschweig , I95I, Z. angew. Math. Mech. 33, 1953, p. 89. 

^■'Htfith the symbols according to W. Bienemann there apply the 
identities Ht = g 0 and. A — K. 
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one obtains after calculation of the integral with the velocity distri- 
butions (16) and (17) the quantity g 0 as a function of K and A. 

According to W. Dienemann, there results 


g 0 (K,A) = g x (A) + Kg 2 (A) 


(27) 


62(A) 

A ^ 1: gq(A) 

g 2 (A) 


= 1a 

2-A 3 + -i-A^ 

15 

140 180 

= 1a 

-iv\2 — 1^/5 1 aJ 

90 

84 560 1^080 A 

_ 5 

3 1 2 1 31 

10 

10 A 15 A 2 140 .4 
^ A 


> 


1 1 _ _1_ 1 _ 
120 A “ 180 a 2 


i_ _i_ . _i_ jl. 
8^° A ^ 


1 _1_ 

180 A 5 

1 _L_ 

3 , 02 k a 5 


(28) 


The function go(A) as a function of A for various values of K is 

represented in figure 5. Table 1 gives a few numerical values of the 
functions gq(A) and g2(A). 


TABLE I.- THE UNIVERSAL FUNCTIONS g 1 (A) AND g 2 (A) 
ACCORDING TO EQUATION (28) 


A 

61(A) 

100g 2 (a) 

0 

0 

0 

.2 

.0053 

.036 

.4 

.0208 

.114 

.6 

.0457 

.205 

• 7 

.0606 

.249 

.8 

.0784 

.291 

• 9 

.0970 

.329 

1.0 

.1175 

.364 

1.2 

.1614 

.423 

1.4 

.2089 

.471 

1.6 

.2589 

.510 

1.8 

.3109 

.541 

2.0 

.3643 

.568 

2.5 

.5021 

.618 

3.0 

.6437 

.651 
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(c) The Momentum Equation for the Meridional Direction 


Further transformation of the momentum equation for the meridional 
direction yields, if one introduces, according to Holstein-Bohlen 1 ^ and 
analogous to ( 23 ) 



v 


= z^ 
dx 


the following differential equation for Z(x) 


(29) 


dZ = F(K,A) 
dx U 


(50) 


Therein 



f 

j— 



F(K,A) = 2' 

f 5 - (2X + Xfo) - Z^ 

^ ^ R dx 

1 + 1 

W f o_ 

> 


exactly as G(K,A) in (26) a universal function of the two parameters 
K and A. Individually, the following relationships apply: 


f 0 (K) = hk - JlL _ JL. _ K 2 
& x 515 94-5 9,072 

(52) 

fi(K) = = J- _ JL. 

s x 10 120 

(55) 

f 2 (K) - 5x * - fl(K) 

(54) 

A ^ 

f 0 (K) K * ~ ■fl 

u XX X 

(55) 

f,(K) = = (h li\ fk £x = ( 2 . K\ f / K x 

5 n u ^ Z y 0 u 5 X r + 6; o (K) 

(56) 

Schlichting, Grenzschicht-Theorie, p. 195 . 
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hQ = 


h. = 22 _ 
6 y 126 


(37) 


The above functions of K are already known from the calculation of the 
boundary layer of the two-dimensional case. ^7 

The connection between Z and K results from (18), with considera- 
tion of (29) and (32), and is 


Kf 0 2 (K) = 


dU 

( v o' 

fu 

dR 

dx + 

u 

/ R 

dx 


(38) 


Taking X = ZU', from (29) , into consideration, one may write this because 
of (32) also in the form 


X 


1 + 


\ U / R U' 



K 

945 


j £_\ 2 

9,072j 


(39) 


In figure 4 the universal functions £ q , ± 2 , f and K are 

represented as functions of X*. At the point of separation, for arbi- 
trary rotational velocity, one will have, because of K = -12, the 
parameter X* = -O.I567. At the stagnation point, without rotation, 

K = 4.716 and X* = O.O 5708 , whereas with rotation the values at the 
stagnation point are dependent on the spin parameter Vq/u (cf . the 

following section). From (38) the form parameter K can be determined 
when Z is given. Furthermore, for the later calculation a connection 
between the parameters A, gQ, 0, Z, and K is needed. There results 
according to (22), ( 23 ), (29)> and (32) as follows 


^0 (K ' A) -/!*oW 


(to) 


The two differential equations (25) for q ( x ) and ( 30 ) for Z(x) 
are two simultaneous differential equations coupled by the universal 
functions G(k,A) and F(K,A). In the case of the nonrotating body, 
v 0 = 0, the coupling is eliminated since then, according to (31).? the 

function F becomes independent of A and remains dependent only on K. 


-^Cf. H. Schlichting, Grenzschicht-Theorie, Chapter XII. 
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In this case, one can first determine Z(x) from (30), and subsequently 
6(x) from ( 25 ). This solution for 0 has - it is true - no physical 
significance. It serves merely for giving the limiting value for 
vanishing speed of rotation. 

(d) The Initial Values at the Stagnation Point 

At the stagnation point where U = 0, the two differential equa- 
tions ( 25 ) and ( 30 ) have a singular value since in both equations on 
the right side the denominator vanishes. In order to obtain at the 
stagnation point initial slopes of finite magnitude, d0/dx and dZ/dx 
finite , the numerators also must disappear in these two equations for 
the stagnation point. This requirement yields the initial values of 
the parameters Kq and Aq at the stagnation point. For the potential 

flow there applies at the stagnation point 


x — >0: U(x) = U 0 'R = aR = 1 

dx 


(41) 


The initial values of the meridional equation are obtained from 
F = 0 according to (31) 


f 30 " 2> b " " *0 


1 + 



= 0 


With 


Voo 2 _ *0 


1 + 



according to ( 38 ), with 



according to ( 36 ), and h^ according to ( 37 ) > there results after a 
brief calculation 
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For a given speed of rotation m/a, this is the first equation between 
the initial values Kq and Aq. For the case without rotation, m = 0, 

the boundary-layer thickness ratio Aq drops out from this equation, 
and an equation for the initial value only remains which reads 

2 + Aiqr + 2 — l q-3 = 0 ' (43) 

105 ^ 2,530 ^ 3,024 ^ ' 

The physically useful solution of this equation is 

i^oj^Q = ^0 = ^Al6 (H) 

as known according to S. Tomotika. 

For the initial values of the azimuthal equation, one obtains from 
G = 0 according to (26) 


2s o( K o ,A o) " cr ( K o^ A o) " 0 

Because of 

& x 2 

0(Kq,Ao) = 4ae 0 = 4ag 00 2 A 0 2_5_ 
according to (24) and 


*0 " 


a^ 2 


+ 



according to (l8) and because of (27) one obtains after a short inter- 
mediate calculation 


Ko 


1 + 



1 

2 [si(Ao) + Ko^^Ad^Ad 2 


(^5) 


For a given speed of rotation m/a, this is the second equation between 
the initial values Kq and Aq. 
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For the case without rotation, to = 0, one obtains from (45) with 
Kq = Kqq = 4.716 according to (44) for the initial value of' 

q = Aqq the equation 


1 - 9.452Aoo" 


' g i(^oo) + ^•7 1 6g 2 (Ado) 


= 0 


(46) 


Hence results with g x (A) and g 2 (A) according to ( 27 ) and ( 28 ) 


Aqq = 0.915 


( 47 ) 


The ratio of the boundary-layer thicknesses A = Sy/& x for the azimuthal 

and meridional velocity distribution therefore lies near 1 which is 
physically plausible. 


The two equations (42) and ( 45 ) now represent, for prescribed 
angular velocity o/a, two equations for the initial values Kq and Aq. 


A solution was obtained by determining from both equations the values of 

v2l 

as a function of Aq for various fixed values Kq. 


1 + 


(8J 


Hence, the initial values indicated in table 2 result. These values are 
presented in figure 5 as a function of cjo/a. It was found that for 
values of o/a > O. 815 , no usable initial values of Kq and Aq exist; 

that is, our method fails for these larger values of oo/a. The limit 
beyond which our calculation method fails coincides with the value K = 12 

-1 O 

of the form parameter^ 0 . The initial values Xq, Xq*, and gQQ deter- 


mined from the initial values Kq and Aq 
and table 2, as a function of o>/a. 


are represented in figure 5 


•L^For K > 12, because of the effect of the centrifugal forces, 
is entirely possible in the present case to obtain velocity profiles 
with u/u > 1. 


it 


1 6 
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TABLE 2.- INITIAL VALUES AT THE STAGNATION POINT 


a 

K o 

A o 

100 Xq 

100 Xq* 

^00 

0 

4.716 

0.915 

5.71 

5.71 

0.0629 

.221 

5 

.908 

5.71 

5-99 

.0632 

.454 

6 

.882 

5.70 

6.89 

.o64o 

.679 

8 

.858 

5.69 

8.32 

.0651 

-785 

10 

.781 

5.69 

9.19 

.0661 

.815 

12 

.726 

5.69 

9-49 

.0664 


Finally we obtain the initial value for Z simply in the following 
manner with Uq' = a 





The initial value for © results with a Q = 2gQ 0 according to (24) as 


% = 


1 800 
2 a 




The expression for the velocity distribution used here (parabola of 
the fourth degree for u and v) is different from that of our former 
calculation^ for the rotating disk in a flow. It must be expected, 
however, that the boundary -layer parameters of the rotating disk in a 
flow should agree approximately with those at the stagnation point of 
the rotating body of revolution if both methods are to yield usable 
results . We give this comparison for the momentum-loss thickness in 
x-direction (8) at the stagnation point and for the meridional component 
of the wall shear stress at the stagnation point. The dimensionless 
momentum-loss thickness at the stagnation point is according to (29) 
with U x _q ' = a 



( 50 ) 


•*■9 gee footnote 5 on page 2. 
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The meridional component of the wall -shear stress at the stagnation 
point (jxoj q = T r 0 is according to (5), (l 6), and (32) 




( 51 ) 


The values calculated accordingly are compared with those of the rotating 
disk in figure 6.20 The agreement up to the validity limit of our cal- 
culation (ca/a = O.815) 'is quite satisfactory. 


Hence we conclude that our present calculation yields satisfactory 
results in the entire range 0 5 : ta/a S O.8I5. 


4 . TORQUE AND FRICTIONAL DRAG 
(a) Torque 


The entire torque of the body of revolution may be easily ascer- 
tained from the results of the boundary-layer calculation in the 
following manner: The contribution of an element of the body of revolution 

with the radius R(x) and the arc length dx is (fig. 7 ) 


3M = -2nR 2 Ty Q dx 


and thus the total torque 



T y 0 E2dx 


( 52 ) 


where signifies the arc length from the stagnation point to the 

point of separation. Taking the momentum theorem for the circumferential 
direction (11 ) into consideration, one obtains 


2 ®Whereas the values for the wall-shear stress could be taken 
directly from the report referred to in footnote 5 (p. 227, table 2), 
the values for the momentum-loss thickness were calculated subsequently 
with application of equation (8) with the velocity distributions indi- 
cated there. 
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M - 2itpo)|R U-9 X yJ ^ 2rtpR^ U^ X y^ 


(53) 


where the subscript A denotes the values at the separation point. 

From the boundary-layer calculation, one knows the value of the momentum 
thickness due to spin at the separation point in the dimensionless form 


^A 

®m 



where is assumed to denote the radius of the maximum cross-sectional 

area. 

If one introduces - in the same manner as for the rotating disk - 
a dimensionless spin coefficient by 


°M ~ 




(55) 


one obtains 



U 

00 



*ktB 


MX 

U ra 


(56) 


where V m = R^ is the circumferential velocity of the maximum cross- 

sectional area. Since, as the completely calculated examples show, the 
dimensionless momentum thickness due to spin B varies at the separa- 
tion point only a little with V m /U M , c M is in first approximation 

proportional to and inversely proportional to the Reynolds 

number j/iy^/v. 

For the case of the rotating disk in a flow, with the radius Rj^ = R, 
one obtains because of R^ = R, = aR from (56) in combination 
with ( 5 ^) 

c m = \ji 



(57) 
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and with the numerical , value 



.. according to (4-9) 





in very good agreement with the former investigation 21 where the numeri- 
cal value is 5 - 17 - 


(h) Frictional Drag 

The frictional drag of the rotating body of revolution may be deter- 
mined by integration of the wall-shear stress components T x 0 ‘ a sur ~ 

face ring element of the body of revolution with the radius R(x) and 
the arc length dx (fig. 7) yields the drag 


dW = SitRT^dS ( 58 ) 


Therein x is the coordinate measured along the body axis. Integration 
from the stagnation point x = 0 to the separation point x^, where 


t xq = °> yields 



Rdx 


(59) 


We shall refer the drag to the maximum cross-sectional area nR^ an d 
define the drag coefficient 


'W 


W 


21cf. footnote 5 on page 2 , equation (4-9a) 


(60) 
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Since we obtain the wall-shear stress in the dimensionless form 


T *o gST 

pUoo 2 V V 


(61) 


we may write for the drag coefficient 




T-ILa^- 

^m ^m 


(62) 


5- EXAMPLES 
(a) Sphere 


As the first example, the friction layer on the rotating sphere was 
calculated. When R m signifies the sphere radius, x the arc length, 

and = cp the center angle measured starting from the stagnation 

point, the radius distribution is 


R(x) = Rm sin cp ■ J ' (65) 

*• 4. 

and the theoretical potential velocity distribution 

U(x) = |u ro sin cp (64) 


The velocity gradient at the stagnation point is 


a = (sn\ = l ik 
V^/Uo 2 R m 

and thus 

ft m\ _ _ 3 cu 

\U J 0 u “ 2 a 


(65) 
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Since, according to the explanations in section 3j> the calculation can 
be carried out only for m/a ^ 0 . 815 , we must limit ourselves to 

VjJu*, ^ 2. 0.815 = 1.22. 

’ 2 / " •»- - 

The solutions are obtained by numerical integration of the two simul- 
taneous differential equations ( 25 ) and ( 50 ) for the two cases = 0 

and 1. The calculation scheme is given in table 3 . The results for 
further values of V m /V m could hence be obtained conveniently by inter- 
polation. The case = 0 (nonrotating sphere) agrees with the case 

of Scholkemeyer^ . The results of the calculation are represented in 
table 4 and figures 8 to 12 . 

TABLE 3.- CALCULATION SCHEME FOR THE SOLUTION OF THE TOO 
SIMULTANEOUS DIFFERENTIAL EQUATIONS (25) AND (30) 


cd prescribed 


■ 

X* 

K 




<p 

X 

R(x) 

dx 

U(x) 

dx 

v 0 = oR 

Z v 

H 

f 0 

' f 2 

f 3 








Initial 

value 

given 

(ect. («3)) 



Fig. 4 
initial 
value 
(table 2) 


Fig. 4 

Fig. k 

<- 



Given 

— body form and potential flow > 

< 


To be calculated 
line by line 


> 



2ff 

S<A 

B 

4 60 


F 

dZ 

dx 

AZ 

z v+l 

G 

de 

dx 


®v+l 

Initial 

value 

(eq. (49)) 

Eq. (24) 

Eq. (40) 

Fig. 5 

■ ( gcA 

A 


Eq. (31) 

Eq. (30) 



Eq. (26) 

Eq. (25) 




^ — — — — ■■ — ■ To be calculated line by line 


TABLE 4 . - POSITION OF SEPARATION POINT AND OF THE TORQUE IN 
DEPENDENCE ON V m U OT FOR THE ROTATING SPHERE IN A FLOW 


Spin 

parameter, 

' v m 
U ~ 

Separation 

point. 

Torque, 

0 

108.2 

9.15 

.25 

108.0 

9-14 

.50 

107.5 

9.06 

.75 

106.2 

9.03 

1.00 

104.9 

8.95 

1.22 

105.5 

8.85 


•Footnote 13 on page 8. 
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Figure 8 gives the variation of the form parameter K of the merid- 
ional component of the velocity distribution in the boundary layer. The 
initial values Kq at the stagnation point are immediately given in 

table 2 with equation (65). At maximum velocity, cp = 90°, K is, 
according to (l8), equal to zero for all V m ^ U m , because in a sphere at 

the point where du/dx = 0, also dR/dx = 0. The value K = -12 gives 
the position of the separation point A. In figure 8 the variation of 
the boundary-layer thickness ratio A = byj S x is also plotted; it always 

lies close to 1 and also changes only little with V^JUoo. Figure 9 shows 
the variation of the momentum thickness due to spin "3xy- The curves for 
various Y m jU x almost coincide. The same is true for the momentum-loss 
thickness ■3 X and the friction-layer thicknesses 8 X and 5y. Figure 10 

shows the variation of the meridional and azimuthal component of the 
wall-shear stress. The meridional component t X q increases with the 

spin coefficient V m /Uco only a little whereas the azimuthal compo- 
nent Tyg in first approximation is proportional to the spin coeffi- 
cient V m M . The position of the separation point as a function of 
the spin coefficient V m ^[J oo is given in table 4. For the nonrotating 
sphere cp A = 108.2°, and for Y m j = 1.22 the separation point shifts 
forward to <p A = 103.5°- This displacement of the separation point 

because of the rotation is due to the effect of the centrifugal forces 
and is, clearly, immediately plausible. For the velocity profiles 
behind the equatorial plane (ep > 90°), the centrifugal forces have the 
effect of an additional pressure increase in flow direction and there- 
fore cause the separation point to shift forward. In figure 11 the 
dimensionless torque coefficient formed according to equation (56) is 
represented as a function of the spin coefficient (Cf. table 4.) 

One sees that the proportionality with is fulfilled with very 

good approximation. Finally, figure 12 shows several velocity profiles 
in photographic reproduction. 

A sphere is rather unsuitable for the comparison of the theoretical 
c al culation with test results, because of the large dead-water zone which 
has the effect that even in the case of the nonrotating sphere the posi- 
tions of the separation point according to theory and to measurement do 
not agree when the boundary-layer calculation is based on the potential- 
theoretical pressure distribution as we have done here. A valid compari- 
son regarding the influence of the rotation on the behavior of the fric- 
tion layer can be made only for a slender body where no noteworthy 
dead-water zone develops. Nevertheless we mention here the measured 
results of S. Luthander and A. Rydberg 2 3. j n figure 13 the drag coef- 
ficient of the sphere in dependence on the Reynolds number Re for 
various values of is given according to these measurements. For 

the nonrotating sphere, V m / = 0, and up to values of to about 

3, the curve c^. against Re shows the characteristic variation with 
23 


Footnote 2 on page 1. 


MACA OM 1415 


23 


the familiar sudden drop at the so-called critical Reynolds number. It 
is known that for Reynolds numbers below the critical Reynolds number 
the friction layer undergoes laminar separation, and for numbers above 
the critical Reynolds number, in contrast,- a turbulent one. In the case 
without rotation, the laminar separation point lies at about cp = 81°, 
the turbulent one, in contrast, at about cp = 110 ° to 120 °. The meas- 
urements with rotation show for Vjn/U^ = 0 to 1.4 a shifting of the 

critical Reynolds numbers toward higher values of Re . This shifting 
of the critical Reynolds number to higher values for small V m /U OT is 

probably brought about by the fact that for V^U^, = 0 the laminar 

separation point is shifted from cp = 8l° to higher cp-values, with the 
separation still remaining laminar, however. Only for higher values of 
Vjn/Uoo, the rotation causes the friction layer to become prematurely tur- 
bulent, and it then has the effect of a trip wire whereby a shifting of 
the critical Reynolds number to lower Reynolds numbers takes place . 

Whereas in our theoretical calculations a forward displacement of 
the separation point occurs, due to the influence of the rotation, the 
measurements for small values of indicate a shifting of the 

separation point toward the rear. On the basis of the effect of the 
centrifugal forces, this must be expected, if one takes into considera- 
tion that in the case without rotation the laminar separation point lies, 
according to theory, behind the equator, according to measurement, how- 
ever, ahead of the equator. In both cases, the separation point is 
shifted toward the equator by the effect of the centrifugal forces as 
is to be expected, at least for small V m ^ !!«,, as long as no premature 
laminar /turbulent transition has been produced by the rotation. 

(b) Bodies With a Base (Half -Bodies ) 

As a second example we shall now treat the so-called half -body 
(body of revolution I) which originates by superposition of a transla- 
tional flow on a three-dimensional source flow. If one denotes by R m 
the largest radius at infinity, the following parametric representation 
for the geometrical data of the body24 is 

-E- = sin ^ (66) 

Rm 2 


^^For these relationships as well as for the numerical calculations 
of section 5a, I am indebted to Dr. E. Truckenbrodt . The example calcu- 
lations of sections yb and c are taken from the thesis of K. H. Gronau, 
1952. 
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£ - tan I / 1 -! sin2 l + *(f'!) (67)25 

Here cp is the angle measured from the forward stagnation point. 

F and E are the incomplete elliptical integrals of the first and 
second kind for the modulus a = 60 °. 

The velocity distribution is 


-U = 2 sin £ ,/l - 2 sin2 £ 
Uoo 2 V 4 2 


( 68 ) 


The form of the body and the velocity distribution are represented in 
figure 14. This figure shows, for various values of the spin param- 
eter the variation of the form parameter K with the distance 

along the body. One sees that already for =» 1.3 only positive 

values of K result . This means that due to the rotation the laminar 
friction layer has become more stable because in the present case the 
centrifugal forces accelerate in the direction of the flow and thus 
have the effect of an additional pressure drop. We shall forego dis- 
cussing here all the results. In figure 15 we have represented the 


torque coefficient 


Vk 




‘Pm 


against the length of the half 


body. Moreover, the asymptotic solution was drawn in for comparison; 
one can derive for it the relationship 


V, 


m 


U„ 




(69) 


Aside from the torque, the frictional drag also was determined. 
Figure l6 presents a compilation of the torque coefficient and of the 
drag coefficient in dependence on the spin parameter for vari- 

ous body lengths Lj 1^- It should be emphasized that the drag coeffi- 
cient is increasing about quadratically with the spin parameter which 
is in qualitative agreement with the test results that have become 
known so far. 


^5e and F signify 



'9/2 


d$ 


fi - k sl 


and E 


sin^ 3 


£ g 

2lb 




sinr 


3 d-d 
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(c) Streamline Bodies 

As further examples we also calculated two streamline bodies of the 
.thickness ratio d/l = 0.2 (bodies of revolution II and III). The body 
shapes and the pertaining velocity distributions were taken from the 

report of A. D. Young and E. Young^ (fig. 17 ). The body of revolution II 
has as a meridional section a normal profile; the body of revolution III, 
in contrast, has a laminar profile with the velocity maximum lying rela- 
tively far downstream. Of the results, figure l 8 shows the torque coef- 
ficient and the frictional drag coefficient as a function of the spin 
parameter V^U^. In both cases, there are not large differences between 

the bodies. For the rest, the variation is s imil ar to that in the case 
of the body with a base. In figure 19, the position of the separation 
points is shown as a function of the spin parameter V m /U 00 . In agreement 

with the values for the rotating sphere (cf. table k), the separation 
point shifts forward with increasing rotational speed. This displacement 
is larger for the body of revolution II than for the body of revolu- 
tion III which is made understandable by the position of the velocity 
maximum. Finally, we gave for the body of revolution II a graphic repre- 
sentation of the velocity distributions in the friction layer for the spin 
parameters = 0 and = 1 (fig. 20). From it one sees that 

ahead of the pressure minimum the meridional velocity component does not 
vary noticeably due to the influence of the rotation whereas between the 
pressure minimum and the separation point the influence of the rotation 
is considerable. 


6. SUMMARY 


A calculation method is given by which the flow about a rotating 
body of revolution in a flow which approaches in the direction of the 
axis of rotation may be determined on the basis of boundary-layer theory. 
The investigations yield a contribution to the aerodynamics of a pro- 
jectile with spin. The calculation is carried out for the laminar 
boundary layer with the aid of the momentum theorem which is stated for 
the meridional and for the circumferential direction. The performance 
of the calculation requires the solution of two ordinary simultaneous 
differential equations of the first order. It yields, in addition to 
the boundary-layer parameters, the frictional drag and the torque as a 
function of the dimensionless spin coefficient V lU = circumferential 

velocity /free -stream velocity. The displacement of the separation point 


26 

A. D. Young, E. Young, "A family of streamline bodies of revolu- 
tion suitable for high-speed and low-drag requirements." ARC Report 220k, 
1951 . 
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with the spin coefficient also is obtained. As examples, the flow about 
a rotating sphere, about a body with a base, and about two streamline 
bodies is treated. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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Figure 2.- Velocity distributions in meridional and in azimuthal direction. 



Figure 3.- The universal function a = &y/s x as a function of gg(A) =5 X y/s. 

according to (27) and (28). 
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Figure 4.- The universal functions f Q = V 8 x > f 2 = f 3 = t x0 / 4 * x /u 

and K as a function of x* according to (32) to (39). 
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Figure 12.- Photographic representation of the velocity distributions 
in the boundary layer of the rotating sphere in a flow, V m /U M = 1, 
A = separation line. 



1.0 2.0 3.0 4.0 5.0X102 


Figure 13.- Drag coefficients of the rotating sphere in a flow as a 
function of the Reynolds number, according to measurements of 
S. Luthander and A. Rydberg. 



Figure 14.- Form parameter K for the rotating body with a base 

flow (body of revolution I). 
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Figure 15.- Torque coefficient c M for the rotating body with a base in 
a flow as a function of the length of the body. Re = U R m /v . 
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Figure 18.- Torque coefficient c ^ and drag coefficient c w for the 

two streamline bodies according to figure 17 as a function of the 
spin parameter V m /U x . Re = U„R m / v . 
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Figure 20.- Velocity distribution in the boundary layer on the body 
of revolution n for the spin parameters V m /U<> = 0 and 
Vm/Uoo = 1.0. Re = U„R m /v. In the lower figure the boundary - 
layer thickness 8 is magnified 1/30/Re -fold. 


